Frictional cooling of positively charged particles 

Daniel Greenwald^^, Allen CaldwelP 

"'Max Planck Institute for Physics, Munich, Germany 



Abstract 

One of the focuses of research and development towards the construction of 
a muon collider is muon beam preparation. Simulation of frictional cool- 
ing shows that it can achieve the desired emittance reduction to produce 
high-luminosity muon beams. We show that for positively charged particles, 
charge exchange interactions necessitate significant changes to schemes pre- 
viously developed for negatively charged particles. We also demonstrate that 
foil-based schemes are not viable for positive particles. 
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1. Introduction 

One of the focuses of research and development towards the construction 
of a muon collider is muon beam creation and preparation. The short lifetime 
of the muon necessitates preparation of a muon beam on time scales shorter 
than a microsecond. To achieve luminosities on the order of 10^^ cm~^ ■ s~^ 
as in the schemes of [ll, 0, 0] , requires a reduction of beam emittance — known 
as beam cooling — by six orders of magnitude. 

Simulations of frictional cooling [^] show that it can achieve the de- 
sired emittance reduction and potentially produce high-luminosity muon 
beams jsf. Previous work on frictional cooling, including both simulation 
and experiment [5|, focused on cooling \i~ beams and assumed the cooling 
mechanism would be identical for beams. However, positively charged 
particles, unlike negatively charged ones, participate in charge exchange pro- 
cesses, and this greatly alters their cooling. In this paper, we present cal- 
culations and simulations of the effects of charge exchange processes on the 
frictional cooling of positive particles. We also show that previous simula- 
tions of foil-based schemes are invalid for positive particles, and rule out the 
viability of such schemes. 
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Figure 1: Stopping power of helium on (velocity scaled from [6|) and the 
acceleration power of an electric field of fixed strength E for a particle of 
constant unit charge as functions of particle energy T. 

2. Frictional Cooling 

Frictional cooling reduces the energy spread and divergence of a beam by 
balancing retarding forces from interactions with a medium, with accelerating 
forces from an electric field to bring the beam to an equilibrium energy. 
Figure [1] shows the stopping power S = (l/p)dT/ds, where dT/ds is the 
energy loss per unit path length, and p is the medium density, for helium 
on (j."^ and the acceleration power (E/p) of an electric field of fixed strength 
E on a particle of constant unit charge. The stopping power is velocity 
scalecH [3] for from the proton data given in [i]. The stopping power on 
\x~ has a similar shape, though it is smaller in magnitude at energies below 
approximately 100 keV [s], 0] . When the accelerating power is larger than 
the stopping power, the particle is accelerated. When the reverse is true, 
the particle is decelerated. At the kinetic energies where the two powers are 



^Except where noted, data for |x+ interactions in this paper are velocity scaled from 
proton data. 
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0.11 
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Table 1: Energy spreads of frictionally cooled beams of ]x~ in H2 gas from 
and yi^ in He; and gradient of the stopping power of H2 on [i" from the 
parameterization in [8| 



equal, particles are at an energy equilibrium. If the stopping power is greater 
than the accelerating power at energies above this equilibrium point and vice 
versa below it, then the point is stable and attractive. 

For a particle of constant charge, two equilibrium points are created: a 
stable one (Teq in figure [1]) at an energy below that at which the stopping 
power peaks, and an unstable one at an energy above (T^q). Particles with 
T < will be brought to Teq, defining the frictional cooling energy region. 

Straggling of energy losses to the medium prevent a beam from becoming 
truly mono-energetic in a frictional cooling scheme, inducing a spread of the 
beam energy around Teq. A study of the frictional cooling of yL~ in found 
that the spread of the energy distribution of a cooled beam is independent 
of the beam's initial spread and, as table [T] shows, decreases with increasing 
gradient of the stopping power. In the table, the gradient is calculated from 
the parameterization of the stopping power for hydrogen on [i~ from j^, 
using the measured parameters given in the paper. 

Such effects are shared by ionization cooling , [m , which operates at 
higher kinetic energies. However, in contrast to ionization cooling, frictional 
cooling, which operates at low kinetic energies, is greatly affected by nuclear 
scattering, which contributes to the spread in kinetic energies of the cooled 
beam to a larger extent than the straggling of energy losses. This is due 
to particles scattering away from the electric field direction, possibly even 
in directions opposed to that of the field; they are then slowed down and 
reaccelerated to the equilibrium energy in the direction of the field. 

3. Low-Energy Stopping Processes 

At high kinetic energies, a projectile slows down in a medium through 
excitation and ionization of the medium atoms. One can neglect the inter- 
actions with the nuclei of the medium and assume that a positively charged 
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projectile is stripped of all its electrons [12|, |13|. Thus the stopping power of 
the medium for positive particles is the same as for negative ones. This is the 
energy region of ionization cooling, and the energy loss here is comparatively 
simpler to model than energy loss at low energies. 

At low energies, projectiles slow down by Coloumbic interactions with 
both the nuclei and the electrons of the medium; furthermore the inter- 
actions of positively charged projectiles involve more than excitation and 
ionization, leading to a difference in the stopping powers for positive and 
negative projectiles. 

Scattering of the projectile particle off of a nucleus results in a loss of 
energy and a change of direction. Though, as described above, nuclear scat- 
tering has the largest impact on the final energy spread of a cooled beam, 
it is not the main mechanism of energy loss. The nuclear stopping power is 
orders of magnitude smaller than the electronic stopping power for all but 
the lowest kinetic energies j^. 

3.1. Electronic Stopping 

The interactions of the projectil^ with the electrons of the stopping 
medium are the dominant mechanisms of energy loss. In the frictional- 
cooling energy region these interactions are the excitation and ionization 
of the medium atoms (X), 

^W+X ^ ^W + X* (la) 
^(+)+X'? ^ ^W + X^+^ + ne" (lb) 

and the capture and loss of an electron by the projectile, 

^(+)+X^ ^ [x^^^+X^+i (Ic) 
^(o)+X^ ^ ^W+X^ + e" (Id) 
^(o)+X^ ^ ^W+X^-\ (le) 

where we have introduced the notation yS'^'' to represent charge states of 
as an ion of muonium (Mu): 

^(+) = Mu+ = , 

^(0) = Mu = ^i+e", and (If) 
yS~^ = Mu" = M.+e-e-. 

It is obvious that stopping of the neutral charge state plays a role in the 
stopping of the projectile. We must also consider processes (fTal) - (fTe|) with 



^The equations to follow are equally valid with |a.+ and Mu replaced by p and H. 
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the replacement of by iJ.^''^ and ij.^*^^ by IJ.*-^-*; and processes involving the 
negative charge state: (llap and fllbp with the replacement of [x^^^ by 
and double electron capture and loss 

We write the cross sections for processes fllap -( IIh|) for \Ji^~^\ ]^^'^\ and ij.*^^-' 
with the notation aqq'(T), denoting the total cross section for the interactions 
taking a muon of charge state q and energy T to a muon of charge state q' 
(accompanied by an energy loss). To be clear, q and q' refer to the charge 
state of the muonium ion (+,0,—), not to the charge of the muon itself, 
which remains positive in these purely electromagnetic interactions. 

The total stopping power of the projectile is the combination of the indi- 
vidual stopping powers of the different charge states, 

S{T) = Y.f'iT)S\T) 

= f^{T)S+{T) + f{T)S%T) + f-{T)S-{T), (2) 

where the are the equilibrium charge state fractions which are the 
solutions to 

^ oc 5: {r'aq,q - /V,,,) = 0, Vg, (3) 
q' 

and 

E/'' = i- (4) 

Since the df^/dx are taken at a fixed T, the are functions of T. For the 
three-state system {m-'-'''^ \Ji^^\ M-*'"''}, the equilibrium charge states are 

r = A^ / Y: A^' , (5a) 

/ g' 

where 

A" = ^ aijajg + n (^iq (5b) 

By the nature of its measurement [sl, the stopping power shown in figure[T] 
is the total stopping power; that is, the left-hand side of equation ([2]). 
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Figure 2: Charge exchange cross sections for in hehum, velocity scaled 
from proton cross sections in 15| resulting in final states with (large 
dashes), (x^'^-' (small dashes), and (dots). The two-state charge-change- 
cycle cross section is shown in solid red. 



4. Effective Charge 

It is clear from section [3?T] that when traveling through matter, [x"*" spends 
some time as Mu and Mu~. This changing of charge state will have a signif- 
icant impact on the frictional cooling process, which involves the restoration 
of energy losses by an electric field according to qE. 

Figure [2] shows the cross sections for the charge exchange processes of yi^ 
in helium obtained by velocity scaling those for protons in [15|. The cross 
sections resulting in yS~'^ charge states are orders of magnitude smaller than 
those taking yS''^ to \\.^^'^ or So ij."*" traveling in helium (and in all the 
materials we will discuss) is nearly in a two-state system {iJ.^"'"^ M-''°^}- This 
simplifies the calculation of the equilibrium charge state fractions to 

/+ = and f = . (6) 

(7+o + c^o+ o-+o + ao+ 

The mean free path for one charge change cycle (fJ-'^^^ — )■ — )■ \Ji^^^) can 
be calculated from the mean free paths for the individual charge exchange 
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Figure 3: Effective charge (line) for in He, H, Ne, and C for tlie tliree- 
state (solid black) and two-state (dashed red) systems. The data points show 
the equilibrium charge state fractions taken from . 



processes, Acy = Ao+ + A+o- This yields a cross section for a charge exchange 
cycle to take place 



a. 



^io^^ (7) 

which we show in figure [2l Since the charge exchange interactions take place 
over distances much shorter than the overall distance the muon travels in 
the frictional cooling scheme, we can approximate the charge of the [x^ by 
an effective charge according to 

r/^E?/"- (8) 



Figure [3] shows the effective charges of [i'^ in helium, hydrogen, and neon 
for both the two-state (772) and three-state (773) systems calculated from 
the charge-exchange cross sections using empirical formulae from fit to 
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measured values [15[. We also show 772 for carbon; 773 could not be calculated 
since the cross sections involving the negative charge state are unknown. 

For comparison to the calculated effective charge, we show experimentally 
measured charge state fractions from [l^] for helium, hydrogen, and neon: 
/"•", which is the same as r]2', and (/"'" — /~), which is the same as 773. The 
calculated effective charge matches very well with the measured data. As 
well, we see that ri2 and rj^ differ only minutely and only over a small range 
of energies. The negative charge state fraction contributes to the three-state 
effective charge at percent level and lower in all three gases. 

It is important to note that helium and neon are the only materials in 
which the effective charge tends to a value of or near unity at low energies. In 
hydrogen and carbon (as well as water, oxygen, and nitrogen, and therefore 
air) the effective charge approaches zero at low energies. 

5. Accelerating Power & Cooling Medium 

The effective charge of a positively charged projectile in the retarding 
medium of a frictional cooling scheme can be absorbed into the accelerat- 
ing power of the electric field. In effect, this makes the accelerating power 
dependent on the projectile's kinetic energy. 

We reformulate here the requirements for frictional cooling stated in sec- 
tion [2j An equilibrium energy must be established by balancing energy loss 
with energy gain, 

5(Teq) = -r/(T,q). (9) 
P 

At energies above Tgq, the stopping power must be greater than the acceler- 
ating power 

5(Teq + e) - -r/(Teq + e) > 0. (10a) 
P 

And at energies below Tgq the accelerating power must be greater than the 
stopping power 

-r/(Teq - e) - 5(T,q - e) > 0. (10b) 
P 

The last two requirements can be combined to one statement about the slope 
of the stopping power relative to that of the accelerating power: 



dS_ _ Edr]_ 
dT rp p dT 

J- cq 

All three requirements can be met when 



= 5'(Teq) - -r/'(Teq) > 0. (lOc) 
T P 

J- CQ 



S'-S^>0. (11) 
ri 
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He H2 Ne C N2 Ar H2O O2 Kr Xe 
without r] 8.4 6.2 14.1 9.6 9.0 7.9 9.0 10.7 9.0 9.0 keV 
with 7] 3.9 0.9 1.1 0.33 0.23 4.5 3.1 keV 



Table 2: Maximum energy at which a stable equilibrium can be established in 
a frictional cooling schemed for with and without accounting for effective 
charge for several stopping media. A blank entry means no equilibrium 
energy can be estabhshed. 

Accounting for the muon's effective charge greatly reduces the maximum 
energy at which the condition of equation f lTT]) is met. Table |2] lists the 
cooling ranges for several materials with and without accounting for effective 
charge. 

The effects of charge exchange processes also reduce the value of T^^, 
the maximum kinetic energy that can be decelerated to Teq. This effect is 
illustrated in figure HJ which shows the frictional-cooling energy region of 
the stopping power of helium on [x"*". Superimposed on the stopping-power 
curve are three electric-field accelerating powers: The black dashed line is 
the naive accelerating power of figure [1], with effective charge neglected, for a 
field strength 560 kVcm^/mg resulting in Teq = 1 keV. The lower red curve 
is the accelerating power for the same electric field strength, accounting for 
effective charge. The equilibrium energy is cut in half, but T^^ remains the 
same. The upper red curve shows the accelerating power accounting for 
effective charge that results in Teq = 1 keV, which requires a field strength 
of 970 kVcm^/mg, and reduces T^^ by a factor of two. 

Achieving the desired equilibrium energy, accounting for effective charge, 
requires a greater electric field strength than is expected in the naive scheme. 
The field strength may in fact have to be an order of magnitude larger, 
depending on what the desired Teq and the stopping medium are (figure [5]). 
On the bright side, the relative slope of the stopping power is generally larger 
when effective charge is accounted for, perhaps causing (according to [3]) the 
final energy spread of a cooled beam to be smaller. As well, the stronger 
electric fields may reduce the spread caused by scattering. 

We performed Monte Carlo simulations of the cooling of (x"*" beams with 
discrete charge exchange interactions, as well as simulations using the effective- 
charge approximation. These conformed to the calculations presented in fig- 
ure O to within a few percent. 

From the simulations we can calculate the spread of the beam energy 
around Teq due solely to charge exchange interactions. Table [1] lists these 
spreads for [i^ in helium at the same equilibrium energies cooled to in Q, 
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Figure 4: Stopping power of helium on iJ."*" (solid) and accelerating power of 
a uniform constant electric field of strength E on ix"*" in helium with (large 
dashes, red) and without (small dashes, black) accounting for effective charge 
for two values of eE (as indicated on the right axis). 

which looked at \x~ in hydrogen. The spread due to charge exchange for 
is smaller than the spreads due to scattering and straggling for \x~ . A 
comparison using in hydrogen is not possible because the maximum Tgq 
in H2 for is 0.9 keV. 

The choice of stopping medium is even further limited than the require- 
ment that the relative slope of the stopping power be positive. This is illus- 
trated by the example of oxygen, in which the relative slope of the stopping 
power is just barely positive over a small region of energies; but it is not 
sufficiently large to clearly establish an equilibrium energy. 

Furthermore, at low energies, oxygen's stopping power is proportional to 
particle velocity, 5* oc T2, but its r] is proportional to T'^, with > |, causing 
the accelerating power to decrease below the stopping power at low energies. 
A particle that experiences a large loss of energy in one interaction — or one 
that scatters into a direction opposite that of the electric field — will not 
reaccelerate to Teq, but rather continue decelerating to thermal energies and 
be lost to the cooling process. This would severely limit the use of such a 
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Teq[keV] 

Figure 5: Ratio of electric field strength required to achieve the desired equi- 
librium energy accounting for effective charge to the field strength required 
without accounting for effective charge. 



medium in a frictional cooling scheme. 

The only viable media for a gaseous frictional cooling scheme for iJ."*" with 
equilibrium energies at or above 1 keV are helium and water vapor. For Teq 
just below 1 keV, hydrogen and neon also become viable; and for Teq of a few 
hundred electron volts, argon and nitrogen are also viable. 



6. Beam Neutralization &: Foils 

The frictional cooling experiment of jsj iterated energy loss and replace- 
ment by means of a series of moderating foils with electric fields between 
them. The experiment and a Monte Carlo simulation in ^ gave promising 
results for such a scheme for the cooling of negative muons. In both publi- 
cations it was posited that the frictional cooling scheme used with could 



also work with iJ."*"; more recently this idea has been revisited in [17[ and [18 
However, the presence of charge exchange interactions greatly limits the yield 
for such a scheme; in fact, leading to a zero yield. 

A foil-based frictional cooling scheme has the benefit of separating energy 
loss, occurring in the foils, from energy restoration, occurring between foils. 
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Since particles are reaccelerated in vacuum, the energy restoration depends 
on the charge state of the particle as it exits the foil. Those exiting in a 
neutral state are blind to the reacceleration field. 

The four studies cited above used and simulated carbon foils in their 
frictional cooling schemes, choosing the reacceleration field strength to pre- 
cisely compensate the energy lost in a foil by a particle at T^q] and assumed 
every particle exits every foil in a charged state. However, for positively 
charged particles, after exiting a foil at mean energy Teq, a beam will be di- 
vided into two populations, a fraction /"'"(Teq) in the ij."*" state and a fraction 
(1 — /"'"(Teq)) in the Mu state. The (x"*" portion of the beam reaccelerates 
to Teq, and upon exiting the next foil in the cooling series, its population 
fraction decreases further to (/"'"(Teq))^. The Mu portion of the beam will 
not have its energy losses to the previous foil restored by the electric field 
and will exit the next foil at an energy below Teq, where the chances of Mu 
atoms losing their electrons are even smaller. Furthermore, Mu atoms that 
have their electrons stripped off, must exit several foils in a row in the pos- 
tive state in order to reaccelerate back to T^q. As the beam passes through 
the array of foils it is neutralized and slowed down. Cooling becomes an 
impractical goal. 

We performed a Monte Carlo simulation of exactlvthis mechanism. The 
foil density was chosen to match the simulations in |j] and the experiment 
in [sl, 5 [xg/cm^ (25-nm-thick foils of 2-g/cm^ amorphous carbon) . Scattering 
and energy-loss straggling were neglected, since they would only worsen the 
beam neutralization described above. Figure [6] shows the yielqj of the foil- 
based scheme as a function of the number of foils traversed for mono-energetic 
beams starting at 9 keV, the maximum Teq possible, and 4 keV, the Teq used 
in After a handful of foils, the yield falls off exponentially with the 
number of foils traversed. After passing through the first foil and energy 
restoration, the mean energy of the beam is greatly reduced from Teq and 
the energy spread is greatly enlarged. The mean energy continues to fall with 
each successive foil. 

We also simulated a scheme in which the foils are coated with frozen 
noble gas to increase the fraction of particles exiting them in the positive 
charge state. This is inspired by the scheme for the production of low-energy 
muons from a surface muon beam at the Paul Scherrer Institute (PSi), which 



uses frozen argon and krypton [19|]. Neither argon nor krypton coatings 



^In the simulation, tracking of a particle is terminated when it does not have sufficient 
kinetic energy to pass completely through a single foil according to Q. The yield is thus 
defined as the fraction of muons surviving passage through the foil. 
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Figure 6: Yield of [x^ in carbon-foil frictional cooling schemes as a function of 
number of foils traversed after reaching Tgq, for both naked foils and frozen- 
neon-coated foils. 



improve the beam neutralization. A neon coating marginally slows beam 
neutralization (shown in figure [H]), but does not prevent the degradation of 
the beam energy. Moreover, we simulated the charge-exchange effects of the 
coatings but neglected their impact on energy loss. This assumes that the 
coatings are nanometers thin; however, the coatings used at PSI are a much 
thicker 200 — 300 nm. Such thicker coatings would hasten beam loss. 



The scheme of [l7j, also based on carbon (graphite) foils, requires the 
beam to pass through order- 100 foils at the equilibrium energy. A beam of 
positively charged muons (or protons) would be unable to survive such a 
scheme. 



7. Conclusion 

Muon collider schemes employing frictional cooling are a viable option 
for collision of multi-TeV lepton beams. Several articles have been published 
with analytical and experimental results for frictional cooling of negatively 
charged particles. Many of these articles have conjectured that the results for 
positively charged particles will be the same as for negatively charged ones. 
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and schemes for the coohng of positive particles have been proposed. A key 
group of physics processes involved in the slowing down of positive particles — 
those changing the charge state — has been neglected in these studies. We 
found that accounting for these processes significantly alters the results for 
positive particles from those for negative ones: The choice of cooling medium 
is greatly limited, such that helium gas becomes the only viable medium; with 
foil-based schemes completely ruled out. The range of equilibrium energies 
for the cooled beam is also greatly limited, with a maximum possible energy 
of approximately 4 keV for jj,"*" (36 keV for protons). And the electric field 
strength required to bring a beam of positive particles to an equilibrium 
energy is significantly greater than the strength required to bring a beam of 
negative ones to the same energy. 
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